Collective Motion of Polarized Dipolar Fermi Gases in the Hydrodynamic Regime 
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Recently, a seminal STIRAP experiment allowed the creation of '"'K^^Rb molecules in the rovi- 
brational ground state [K.-K. Ni et al., Science 322, 231 (2008)]. In order to describe such a 
polarized dipolar Fermi gas in the hydrodynamic regime, we work out a variational time-dependent 
Hartree-Fock approach. With this we calculate dynamical properties of such a system as, for in- 
stance, the frequencies of the low-lying excitations and the time-of-flight expansion. We find that 
the dipole-dipole interaction induces anisotropic breathing oscillations in momentum space. In ad- 
dition, after release from the trap, the momentum distribution becomes asymptotically isotropic, 
while the particle density becomes anisotropic. 

PACS numbers: 21.60. Jz,67.85.Lm 



Even before the realization of Bose-Einstein condensa- 
tion (EEC) with ^^Cr [l|, much experimental and theo- 
retical interest has been dedicated to ultracold quantum 
gases interacting through the long-range and anisotropic 
dipole-dipole interaction (DDI) For bosonic dipolar 
particles, the starting point of the theoretical investiga- 
tions was the construction of a corresponding pseudo- 
potential by Yi and You [3|. After that, an exact so- 
lution of the Gross-Pitaevskii equation in the Thomas- 
Fermi regime was found for cylinder-symmetric traps 
Moreover, the DDI has been shown to shift the 
BEC critical teniperature in a characteristic way in po- 
larized systems Q and to give rise to the Einstein-de- 
Haas effect, when spinorial degrees of freedom are con- 
sidered . From the experimental point of view, time- 
of-flight (TOF) techniques demonstrated both the first 
DDI-signature through small mechanical effects (7^ as 
well as strong dipolar effects in quantum ferrofluids Q. 
Furthermore, the shape of the trap was manipulated to 
stabilize a purely dipolar BEC against collapse [9,] . 
Concerning fermionic dipolar systems, recent theoreti- 
cal studies have considered interesting properties of ho- 
mogeneous gases such as zero sound [Hi], Berezinskii- 
Kosterlitz-Thoules phase transition and nematic 

phases [l^. In harmonically trapped systenis, amaz- 
ing predictions like anisotropic superfluidity [13|, frac- 
tional quantum Hall physics and Wigner crystal- 
lization [1^ have been made. With respect to experi- 
mental investigations, the most promising atomic candi- 
date is the fermionic chromium isotope ^^Cr [l^, which 
has a magnetic moment of m = 6 Bohr magnetons. For 
these atoms, calculations of equilibrium properties have 
shown that the DDI is only a small perturbation [13, [3 ■ 
However, by applying a stimulated Raman adiabatic pas- 



* Electronic address: ' lima@physik.fu-berlin.de' 
t Electronic address: , axel. pelster@fu-berlin.de. 



sage (STIRAP) process, it has recently been achieved to 
cool and trap ^°K*^Rb molecules into their rovibrational 
ground-state, where they possess an electric dipole mo- 
ment oi d = 0.566 Debye 19-22]. Due to the resulting 
strong DDI a considerable deformation of the momen- 
tum distribution is expected JJi jl8j] . Once these systems 
would have been further cooled into the quantum degen- 
erate regime, the main task will be to identify unam- 
biguously the presence of the DDI. In this respect, TOF 
experiments and oscillation frequency measurements rep- 
resent the most fundamental diagnostic tools in the field 
of ultracold quantum gases. Their outcomes reveal im- 
portant information on the nature of the system under in- 
vestigation. They differ drastically depending on whether 
the system is in the collisionless (CL) regime, where colli- 
sion rates are small, or in the hydrodynamic (HD) regime, 
where collisions take place so often that they lead to local 
equilibrium. To date, investigations of dynamical prop- 
erties of trapped dipolar Fermi gases have either been 
restricted to the CL regime [2^ or excluded a deforma- 
tion of the momentum distribution in the HD regime [13] ■ 
Since the experiments with ultracold polar molecules are 
performed under strong dipolar interactions, one should 
expect them to lead the system into the HD regime, and 
thus an analysis allowing for an anisotropy in the mo- 
mentum distribution has to be carried out. In this letter, 
we shall use a variational time-dependent Hartree-Fock 
approach to address this question. 

Consider N spin-polarized fermionic dipoles of mass 
M trapped in a cylinder-symmetric harmonic potential 
C/tr(x) = Muj'^ (x^ + A^z^) /2 with trap anisotropy 
A at ultralow temperatures. Since the Pauli principle 
inhibits a contact interaction, they interact dominantly 
through DDI. As we assume that the fermionic cloud is 
polarized along the symmetry axis of the trap, the DDI 



potential reads Vdd(x) 
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For mag- 



netic dipole moments m the DDI is characterized by 



Cdd — 



whereas for electric moments d we have 
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C'dd = iTTcP . In the following we restrict ourselves to the 
normal phase in the limit T ^ because the critical tem- 
perature for superfluidity is very low, depending expo- 
nentially on add — MCdd/ {4:71^,"^) Furthermore, this 
limit is restricted by the HD requirement that the relax- 
ation time TR is small in comparison with the time scale 
1/w defined by the average trap frequency aJ = {aj'^WzY^^- 
The necessity of a HD approach can be inferred as fol- 
lows. As Tn is not kown for dipolar interactions, we esti- 
mate it by assuming the DDI to be equivalent to a con- 
tact interaction with scattering length a^d- Then we use 
the fact that for a two-component, degenerate, normal 
Fermi gas with contact interaction one has (lutji)~^ — 
{N'^/^aM^/Muj/hfF{TlTY), where F{T/Tf) is of the 
order 0.1 in the quantum temperature regime (see, e.g., 
[25j). Thus, we expect for the one-component, dipolar 
gas to enter the HD regime for A^^/^e^d ^ 1, with the 
dimensionless parameter Cdd = Cdd{M^uj/h^)2 /An 
measuring the strength of the DDI. In the current set- 
up of Ref. nil one has 4 x 10^ '^^K^'^Rh molecules with 
radial trapping frequency of lu^ = tOy ~ 2tt x 175 Hz. 
Assuming an average trap frequency of that value yields 
at least edd « 5.3 and {uJtr)-^ w 0.1 x (A^^/^edd)^ ~ 96, 
which drives the system into the HD regime. 

In this letter we work out a time-dependent 
Hartrec-Fock approach by extremizing the action 
A = J dt{-^\ih-§^ - where ,XN,t) = 

{xi,--- ,xn\'^) is a Slater determinant and H denotes 
the underlying Hamilton operator. In order to de- 
scribe the HD regime, we follow a standard proce- 
dure of nuclear physics and assume that frequent 
particle collisions assure that all one-particle orbitals 
have the same local phase x(2:,t), yielding the ve- 
locity field V = Vx. Thus, we can factorize out 
the phases and define a Slater determinant through 
*o(xi,--- ,a:jv,i) = e-^^^^.x(--*)/''^'(a;i,--- ,xjv,i), 
which contains only the moduli of the one-particle or- 
bitals and, therefore, is invariant under time rever- 
sal. This yields a time-even one-body density matrix 
Poix,x';t) = e-**4x(-.*)-x(^',t)]/K^(a.^a.'.i) [27]. With 
this the action reduces to 
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where po(x;t) = po{x,x;t) denotes the particle density 
and (5'o|-ff|5'o) consists of the kinetic energy i?ki, the 
trapping potential i?tr, and the interaction. The latter 
is divided into the direct or Hartree term E^^ and the 
exchange or Fock term E^^. Due to the exchange term, 

the ground-state energy (^fol-ff l^'o) is not a function of 
the particle density po{x;t) alone, but also contains the 
non-diagonal part po{x, x';t). 

As it is not possible to solve analytically the result- 
ing Euler-Lagrange equations for x(a;,t) and po{x,x';t), 
we propose here a variational extremization of the action. 
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FIG. 1: (Color Online) Spatial aspect ratio for different trap 
anisotropies A; the upper (continuous) branches correspond 
to a local minimum of the mean-field energy and the lower 
(dotted) branches to a maximum. Notice that the value of 
Edd in which two branches meet, i.e., tdd') decreases slower for 
lower values of A. The vertical line marks the estimated value 
of the interaction strength for ''"K^^'^Rb molecules edd ~ 5.3. 



To this end, we express each energy contribution in terms 
of the Wigner transform of the one-body density matrix 
i'o(X,k;t) = /d^spo (X + |,X- §;i) e-*'^•^ The ki- 
netic and trapping energy arc then given by 



E, 



ki/tr 



^iyo(x,k;t)eki/tr(x,k) (2) 
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with Eki = /i^k^/2M and etr = t/ti(x), respectively. The 
direct term, which accounts for the deformation of the 
particle density, and the exchange term, which is related 
to the momentum space deformation, read 

fd^xd^kd^x'd^k' ,,,,,, n , , , , n 
^dd - / KTTT^a !^o(x,k;t)ydd(x-x')!^o(x',k';0, 
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i.o(X,k;t)Kid(s)j^o(X,k';i) 
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At this point, we adopt the variational ansatz 
x(a;,t) — [ax{t){x^ + y"^) + az{t)z'^] /2 for the phase and 

Vq (^X, Jt, tj - U - -j^^ - -j^jj^ - -^jj^ - -j^^ j 

for the Wigner phase space function with G(-) being the 
step function. We are now in the position to extremize 
the action ([TJ with respect to the time-dependent vari- 
ational parameters ai{t) for the phase as well as Ri{t) 
and Ki{t) for the Thomas-Fermi radii and the Fermi mo- 
menta. At first, one obtains = Ri/Ri, which is used 
to eliminate the parameters ai from the rest of the for- 
malism. Under conservation of the particle number 

RlR^klk, - 1, (4) 

the equations of motion for the Thomas-Fermi radii read 

1 d^R., -j^. + + edd A( j^. , j^. , i^.) , (5) 
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FIG. 2: (Color Online) The lower (red) curve shows the ratio 
of the amplitudes (x/Cz Eis a function of edd for A = 5. For 
comparison, the equilibrium aspect ratio in momentum space 
against edd for A = 5 is depicted by the upper (blue) curve. 



Here we use • to represent the quantity • expressed in 
units of the non-interacting Thomas- Fermi radius Rf''^ — 
^2Ef / Mijj'j and the Fermi momentum Kp = \f2WpJW 

1 /3 

with the Fermi energy Ef = (67V) '-"fi-w. The auxihary 
functions are defined according to 
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with the numerical constant Cd = 323/6.5 7 ^2 
Furthermore, the anisotropy function 



0.2791. 



fs{x) = < l-^"" (l_2,2)3/2 , X ^ i 

I 0; X = 1 



(7) 



decreases monotonically from 1 at x = to —2 at a; = 00, 
passing through zero at x = 1 d, Q. In addition, the 
variational parameters are restricted to obey 



eddC 



Rz,K,,K, 



with C 



3cd 



1 



(8) 



This equa- 



tion can be traced back to the exchange term and shows 
explicitly that a non-zero edd implies a deformed momen- 
tum distribution K^^^ for finite Rx , Rz as was first 
pointed out in Ref. [17|. 

Equations dS]) govern the static as well as dy- 

namic properties of a polarized dipolar Fermi gas in the 
HD regime and represent the main result of this letter. 
They determine the temporal evolution of both the spa- 
tial and the momentum distribution of a dipolar Fermi 
gas which are directly experimentally accessible via TOF 
techniques. The static solutions agree precisely with the 
ones obtained before in Refs. [13, U^- In Fig. [T] we 
present our findings for the spatial aspect ratio as a 



FIG. 3; (Color Online) Excitation frequencies for A = 5 as 
functions of the DDI-strength edd. The upper blue (lower 
red) curve represents the monopole (quadrupole) frequency 
r2+ (f^-). The dashed (dotted) horizontal line represents the 
monopole (quadrupole) frequency of the non-interacting gas 
from Ref. [2y|. Inset: Mono- (blue) and quadrupole (red) 
oscillation frequencies of the dipolar Fermi gas normalized 
by the non-interacting values from Ref. [2^ against the trap 
aspect ratio A for different values of the dipolar strength edd- 
The dashed (solid) curves are for edd = 0.5 (edd ~ 1.0). 



function of the dipolar strength edd- The characteris- 
tic feature is that a minimal value of A is required for 
stabilizing a system with a given edd- Thus, for future 
experiments with ^"^K^^Rb molecules in the quantum de- 
generate regime one should choose the anisotropy A to 
be larger than the minimal value Amin ~ 3.89 in order 
to render the system stable against collapse. Amazingly, 
the minimum value of A supporting stabihty, decreases 
slowly and samples with A ~ 0.05 are stable if edd ~ 1-6- 

Having summarized the most important aspects of the 
static solutions, we turn now to their dynamical prop- 
erties. In a cylinder-symmetric system the mono- and 
quadrupole low-lying oscillation modes couple to each 
other. In order to obtain the frequency of these modes 
in the HD regime, we expand the radii and momenta 
around their respective equilibrium values according to 
R, = R,{0) + me'^\K, = K,{0) + Qe'^\ where r?, (Q) de- 
notes a small oscillation amplitude in the i-th direction in 
real (momentum) space and Q represents the oscillation 
frequency. Inserting these into the equations of motion 
®, a linearization yields at first for the ratio of 
the momentum amplitudes 



^ Kx ki + ki - e^^kzdc/dkz 
C. kz 2ki - EMkzdc/dkz 



(9) 



where all terms are evaluated at equilibrium. This quan- 
tity is plotted against edd for A = 5 in the red (lower) 
curve in Fig. [2] and is compared to the corresponding 
equilibrium momentum aspect ratio (blue, upper curve). 
Setting C = 0, i.e., removing the exchange term, one 
has — Cz, whereas for non-zero C, the ratio Cx/Cz de- 
creases monotonically from 1 to about 0.28 in the inter- 
val < edd < fidd' ~ 7.34. This shows that the exchange 
term induces characteristic anisotropic breathing oscilla- 
tions in momentum space, which can be regarded as a 
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FIG. 4: (Color Online) Cloud aspect ratio in TOF expansion 
for A = 5 with tdd = 1,3 and 7 (continuous, top to bottom). 
The dashed curves depict the corresponding momentum as- 
pect ratios. Inset: asymptotic cloud aspect ratio against edd- 

trademark sign of the DDI in fermionic quantum gases. 
Eliminating the momentum amplitudes Q yields a re- 
duced linear homogeneous system for the spatial ampli- 
tudes rji. Demanding non-trivial solutions yields an ex- 
plicit but lengthy result for the monopole (quadrupole) 
oscillation frequency (51-) which depends via the 
equilibrium values of the Thomas-Fermi radii and the 
Fermi momenta upon the trap anharmonicity A and the 
dipolar strength edd- In the special case of an ideal 
Fermi gas, i.e. edd — 0, the oscillation frequencies V,± 

reduce to the correct non-interacting values fl^ = 

(5 + 4A2 ± ^25 - 32A2 -|- IGA^) /3, which were first 

obtained for A = 1 in Ref. [H] and for A 1 in Ref. 
Fig. [3] shows the oscillation frequencies of the mono- 
(blue) and quadrupole (red) modes plotted against edd 
for A = 5. As edd becomes larger, we find that the 
monopole frequency increases and that the quadrupole 
frequency decreases, vanishing at e^^* « 7.34, the same 
value for which the system becomes unstable (see Fig.[T|). 
The inset of Fig. [3] shows how the frequencies depend on 
the anisotropy A for edd — 0.5 (dashed) and edd = 1-0 



(continuous). It turns out that the quadrupole frequen- 
cies are larger than in the non-interacting case for A < 1 
and smaller for A > 1, while the contrary is true for 
the monopole modes. This behaviour agrees qualitatively 
with dipolar BECs 4]. 

It remains to study the TOF expansion of a dipo- 
lar Fermi gas. This is done by numerically solving the 
Eqs. (H])-®, ([S]), while removing the trap frequencies. 
The results are presented in Fig. |4l where the spatial 
and momentum aspect ratios are plotted as functions of 
time in units of uj~ ^ at A = 5 for different edd • The char- 
acteristic of the hydrodynamic regime is that the asymp- 
totic value of the aspect ratio in real space depends on 
edd, while local equilibrium renders the momentum dis- 
tribution asymptotically isotropic. We can estimate the 
validity of these results if we assume the previous HD 
criterion to be valid also during the expansion. Since the 
equations of motion imply d?Ri/dt^ = for large times, 
yielding [RlR^fl^ - aJt, one obtains a HD expansion 
provided (wi)^ • TUtb, ^ 1. For KRb molecules, the ex- 
pansion is HD only for small times uJt 10, whereas for 
molecules like LiCs with d ~ 5.5 Debye, the expansion is 
HD for ZJt < 1000. 

In the present letter we have investigated both low- 
lying oscillation frequencies and TOF expansion data for 
a polarized dipolar Fermi gas through a hydrodynamic 
approach. Our findings have revealed different finger- 
prints of a strong DDI. We have estimated the validity 
of our results and found strong evidence for hydrody- 
namic behavior also in the absence of superfluidity. The 
prospects for observing normal dipolar hydrodynamics 
in the quantum degenerate regime are enhanced by tight 
traps and the recently obtained large dipole moments. 
We acknowledge support from the DA AD, the Innova- 
tionsfond FU-Berlin, and from the DFG in SFB/TR12. 
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